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The Atmosphere Revitalization Recovery and Environmental Monitoring (ARREM) 
project was initiated in September of 2011 as part of the Advanced Exploration Systems 
(AES) program. Under the ARREM project, testing of sub -scale and full-scale systems has 
been combined with multiphysics computer simulations for evaluation and optimization of 
subsystem approaches. In particular, this paper describes the testing and modeling of 
various subsystems of the carbon dioxide removal assembly (CDRA). The goal is a full 
system predictive model of CDRA to guide system optimization and development. The 
development of the C0 2 removal and associated air-drying subsystem hardware under the 
ARREM project is discussed in a companion paper. 




Nomenclature 

a 0 

= 

constant in single-species Toth equation, mol kg' 1 kPa' 1 

E 

= 

constant in single-species Toth equation, K 

c 0 

= 

constant in single-species Toth equation, K 

bo 

= 

constant in single-species Toth equation, kPa' 1 

to 

= 

constant in single-species Toth equation 

a 

= 

single-species Toth equation parameter, mol kg' 1 kPa' 1 

b 

= 

single-species Toth equation parameter, kPa' 1 

h 

= 

single-species Toth equation parameter 

B 

= 

length of test bed filled with sorbent, m 

Rb 

= 

radius of sorbet-filled test bed cylinder, m 

d 

= 

thickness of test bed cylinder, m 

^sorbent 

= 

porosity in the sorbent bed 

P 

= 

density, kg m' 3 

D pellets 

= 

mean diameter of sorbent pellets, m 

^sorbent 

= 

permeability through the sorbent beds, m 2 

Kf 

= 

permeability correction factor 

U 

= 

velocity, m s' 1 

Min 

= 

flow velocity at inlet, m s' 1 

P 

= 

pressure, Pa 

R s 

= 

specific gas constant, J kg' 1 K 1 

Rg 

= 

universal gas constant, J mol' 1 K 1 

M-mix 

= 

molecular mass of gas mixture, kg mol' 1 

Mu 

= 

viscosity, Pa-s 

T 

= 

temperature, K 

T gin (0 

= 

measured time-dependent inlet temperature of gas, K 

Rvapin(0 

= 

measured time-dependent inlet sorbate partial pressure, Pa 

P 

1 vap 

= 

sorbate partial pressure, Pa 

Pout 

= 

calculated absolute pressure at bed exit, Pa 

Pin 

= 

measured absolute pressure at bed inlet, Pa 

AP 

= 

measured pressure drop across bed, Pa 

F 

= 

flowrate, SLPM 

A 

= 

ratio of total sorbent surface area to sorbent volume, m' 1 

A f 

= 

free flow area of can inlet, m 2 

A c 

= 

cross sectional area of the can, m 2 
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Po 

To 

Po 

Pi 

c 

q 

q* 

Gin 

D x 

x 

t 

km 

C P 

k 

h 

dU 

f 

Dab 

Sc 

Re 

Pe 

Nu 

Pr 

ARREM 

AES 

SG 

LPM 

SLPM 

RHS 

CBT 


reference pressure, Pa 

reference temperature, K 

outer perimeter of test can, m 

inner perimeter of test can, m 

concentration of sorbate, mol m' 3 

pellet loading, mol m' 3 

equilibrium pellet loading, mol m' 3 

concentration of sorbate at inlet, mol m" 3 

axial dispersion coefficient, m 2 s' 1 

spatial coordinate, m 

temporal coordinate, s 

mass transfer coefficient, s" 1 

heat capacity, J kg' 1 K 1 

thermal conductivity, W m 1 K 1 

heat transfer coefficient, W m' 2 K 1 

differential heat of adsorption, kJ mol' 1 

directionally dependent coefficient for Pe 

binary mass diffusion coefficient, m 2 s' 1 

Schmidt number 

Reynolds number 

Peclet number 

Nusselt number 

Prandtl number 

Atmosphere Revitalization Recovery and Environmental Monitoring 
Advanced Exploration Systems 
Silica Gel 
liters per minute 

standard (1 atm, 0 °C) liters per minute 

right hand side 

cylindrical breakthrough test 


Subscripts: 

x 

r 

s 

g 

A 

c 


component in the x (axial) direction 
component in the r (radial) direction 
pertaining to the sorbent 
pertaining to the gas 
pertaining to the ambient environment 
pertaining to the can housing 


I. Introduction 

P REDICTIVE simulation tools are being developed to reduce the hardware testing requirements of the 
Atmosphere Revitalization Recovery and Environmental Monitoring (ARREM) project as part of NASA’s 
Advanced Exploration Systems (AES) program 1,2 . Although sub-scale testing is required to establish the predictive 
capability of the simulations, the much greater cost of extensive full-scale testing can be limited to that required for 
the confirmation of analytical design optimization studies. Once predicative capability is established, geometric 
reconfiguration of a model is usually straightforward. A predictive simulation capability provides improved 
understanding of complex processes since process conditions (temperature, pressure, concentrations, etc.) may be 
examined anywhere in the sorption column. Weaknesses in a prototype design can be readily identified and 
improvements tested via simulation. Finally, the predictive simulation provides a powerful tool for virtual 
troubleshooting of deployed flight hardware. Here, we discuss using the COMSOL Multiphysics code 3 to model in 
detail - and predictively - experiments that are similar to the desiccant subcomponent of a full Carbon Dioxide 
Removal Assembly (CDRA) system. 
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II. Design Approach Using COMSOL 

Adsorption in packed fixed beds of pelletized sorbents is presently the primary means of gas separation for 
atmosphere revitalization systems. However, structured sorbents are emerging as a new approach to sorbent 
systems. Structured sorbents are produced as monoliths, with an open structure for airflow, or by fixing sorbents on 
an inert substrate such as paper-like honeycomb structures or expanded metal sheets. A well designed structured 
sorbent is not as subject to attrition (e.g., due to fines or dust generation) as a packed bed. Also, by using a 



Figure 1. CBT Test Bed. 


thermally conductive substrate, the heat of adsorption 
can be transferred out of the bed, possibly to the cold 
desorbing bed if geometry permits. However, 
structured sorbents must be evaluated to determine 
their applicability to commercial processes and space 
flight. It must be shown that, in addition to providing a 
more robust solution, the resource requirements (i.e., 
weight, power, volume, etc.) are similar to, if not less 
than, the state-of-the-art packed bed configuration. 

An accurate assessment of structured sorbents and 
comparison with packed bed designs is desirable; 
experimental results so far show unanticipated 
variation in packed bed breakthrough for identical beds 
held under the same conditions. It is suspected that 
small packing irregularities can propagate downstream 
in large beds and impact process efficiency. This 
indicates a margin of error inherent in packed bed 
fabrication and thus a likely superiority of structured 
sorbents for process efficiency and control. This paper 
discusses fully predictive modeling results using 
COMSOL’ s Multiphysics code for a geometrically 
simple fixed bed design. Insights learned from this 
work will be used in future modeling of the entire 
CDRA system. 

For the bulk separation of C0 2 and H 2 0, 
temperature changes due to the heat of adsorption are 
significant, requiring the simulation of the heat balance 
equations through the beds and the housing, as well as 
the equations for sorption processes and fluid flow. For 
columns with small tube diameter to pellet diameter 
ratios, as encountered in internally heated columns, 
flow channeling along the column wall can have a 
strong influence on overall performance. In non- 
cylindrical flow, the influence is great enough to 
necessitate the use of 3-D simulations. Here, with over 
a dozen pellets per cylinder diameter, 1-D models 
should prove accurate enough for predictively driven 
system design. 


III. Cylindrical Breakthrough Test 


A. Description 

The CDRA requires a water-saving bulk drying stage prior to downstream C0 2 removal. The primary goal is to 
continuously remove at least 80% and up to 100% of water vapor from a process air stream. The Cylindrical 
Breakthrough Test (CBT) was constructed to compare sorption kinetics for various sorbent and sorbate pairs. The 
tests consist of flowing a constant amount of sorbate and carrier gas through a fixed bed containing a regenerated (or 
dried out) sorbent. After some period of time (the ‘breakthrough time’), the sorbate is detected at the bed exit. A 
plot of the sorbate concentration or partial pressure versus time is the breakthrough curve. Axially routed 
thermocouples are used to acquire temperature curves inside and outside of the bed as well as before and after the 
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packed bed. The three thermocouples inside the bed are on axis and are located 1 inch inside the inlet and exit and 
in the axial middle of the sorbent-filled bed. Absolute and differential pressure is measured at the column inlet and 
across the column, respectively. 

The test apparatus, shown in Fig. 1, was designed to have relatively low mass (to reduce regeneration time) and 
good axial symmetry (to reduce system complexity). The test bed, packed with regenerated sorbent pellets, is 
insulated to minimize system heat loss. Mass flow controllers are used to blend N 2 , the carrier gas, with the desired 
partial pressure of C0 2 or H 2 0. For C0 2 tests, Sable Systems C0 2 analyzers provide inlet and outlet C0 2 partial 
pressure readings, while for H 2 0 tests, a Sable Systems Dew Point Generator provides humidity control and 
Edgetech Dewmaster dew point analyzers provide inlet and outlet dew point measurements. 

A schematic of the entire test setup is shown in Fig. 2. The cylindrical column of sorbent has a bed length, B , of 
6.5 inches and a radius, R B , of 1.74 cm. The mass of regenerated sorbent in the test bed is measured so as to 
determine the mean porosity of the bed. The A1 6061 housing is d=0.065 inches thick and extends for 6 inches 
upstream and downstream of the sorbent. The sorbent is held in place inside the housing using spring-loaded plates 
and fine mesh screens. Here, we focus on the CBT adsorption tests, where the sorbent starts fully regenerated. 
Experimental data from the CBT will be used to validate the simulation process so that simulation-driven 
optimization may be used alongside conventional design methods to perfect CDRA sub-component design and 
testing. 

B. Models 

1-D models of the CBT were constructed using the COMSOL Multiphysics code using Domain ODEs and 
DAEs. Only the sorbent-containing part of the bed is modeled. A constant porosity, e sorhent , was used, such that, 
together with the known constant density of the sorbent, p s , the proper measured total sorbent mass is recovered. 
Using the measured mean particle diameter of the sorbent pellets, D pe u ets , the local permeability within each cell is 
then found 4 : 
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K sorbent 


K f e sorbent * ^pellets 

5d5 


(i). 


However, the correct diameter to use in Eq.(l) is a volume -weighted averaged such as the Sauter mean, but this is 
not as well known as D pellets for all sorbents, but it is generally larger than D peaets . The factor, Kf, is adjusted to 
reproduce the measured pressure drop through the test bed for a given test (so it is not a free parameter). Since 
Eq.(l) was developed for typical sandstone packing, some sorbents that are very spherical and smooth have rather 
large values of Kf. Note that for the CBT, the measured pressure drop is close to the analytical pressure drop 
predicted by the Ergun equation. Flow PDEs are derived based on simplification of the Navier-Stokes and 
Brinkman equations while assuming compressible flow; these are combined with PDEs for the loading, based on the 
Toth equations, the sorbate concentration, and the temperature of the gas mixture, sorbent pellets, and housing walls. 

The first PDE finds the 1-D interstitial velocity, u : 

p du = _fdP_ + m \ _ u (2) 

R s T g dt Vdx KsorbenJ P s dt 

where P is the gas pressure, R s is the mixture’s specific gas constant (=R g /M mix ), fi is the gas viscosity, and T g is the 
temperature of the flowing gas mixture. Eq.(2) is essentially a 1 st order Ergun 5 equation that reduces to Darcy’s Law 
in steady state. The 2 nd order Ergun term, somewhat akin to a Forchheimer drag term proportional to the square of 
the fluid velocity, was initially added to the RHS of Eq.(2), but the impact was negligible and so removed for the 
work discussed here. Here, for simplicity the molecular mass of the mixture, M^, is assumed to be that of the pure 
carrier gas, N 2 . Note that the velocity, which is in the axial direction, is not compensated for loss of sorbate since 
the sorbate gas mole fraction is «1. In COMSOL, the entirety of the RHS of Eq.(2) is written as a source term in 
the General Form PDE; the pressure gradient cannot be written as a conservative flux term or the boundary 
conditions become over-constrained. Also, Eq.(2) assumes the ideal gas law for the density of the gas mixture, so 
that: 


^ R S Tg (3) ’ 


The second PDE solves for the gas pressure, P, using the continuity equation and Eq.(3): 

f sorbent dP d ( UP \ 6 sorbent p dT p _ 

RJ g ~dt d^\RJ r g )~ R s T g 2 ~dF~ 1 

In COMSOL, in the Coefficient Form PDE, the 1 st term of Eq.(4) is a Mass Coefficient term, the 2 nd is a Convection 
term, and the 3 rd is an Absorption Coefficient term. The outlet pressure boundary condition is of the Dirichlet type, 
using the measured outlet pressure from a given test, P out = Pin-AP. The inlet uses a mass flux boundary condition 
based on the measured constant standard flow rate, F, for a given test: 


/m — 


Pq F 

P-s'Pq € sorbent A f 


( 5 ), 


where Af= 7iR B 2 is the inlet free flow area and P 0 and T 0 are the reference pressure and temperature, respectively, at 
which the flowrate is defined. Thus, the 1 st term in Eq.(5) is a reference density. Here, P 0 = 1 atm and T 0 = 0°C. 
Since this is a 1-D ‘plug flow’ model, there is no gradient of pressure, velocity, concentration, temperature, or 
porosity in the radial direction. 

The third and fourth PDEs are coupled and solve for sorbate concentration, c, and pellet loading, q , respectively. 
The latter applies the Linear Driving Force model 6 . Together these two PDEs are referred to as the ‘Mass Balance’ 
equations. The General Form PDEs are: 


dc (1 6 sorbent ) dq d 

dt € sorbent dt dx 


( dc\ 
\ Dx fa) ~ 


dx 


( uc ) ( 6 ) 
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— = k m (q, - q ) (7), 

where D x is the local time-dependent axial mass dispersion coefficient, q* is the equilibrium loading from the Toth 
isotherms, and k m is the constant mass transfer coefficient (see below). For the concentration, a zero gradient 
Constraint: 


dc 


dx 


X=B 


= 0 ( 8 ), 


is used at the exit, while a molar volume flux boundary condition is used at the inlet: 

= ( u in c in - uc ) (9), 

x=0 



where u in = is the inlet interstitial velocity and c in = P V apin(t)/Rg/Tgm(0> is the sorbate 

concentration at the upstream inlet to the bed. P va pin(0 and T gin (t) are the measured time -dependent sorbate partial 
pressure and gas temperature, respectively, for each experiment. Since T gin (t) is measured far upstream of the start 
of the sorbent in the bed, the COMSOL inlet value is the measured value minus an offset of 1.2°C*8SLPM/F due to 
losses along the uninsulated piping. That is, this ensures that T gin (£=0) actually equals the T(£=0,v=0) inside the bed 
for any given experiment. There are no explicit spatial boundary conditions required for Eq.(7) since there is no 
loading outside of the sorbent bed. In COMSOL, it is necessary to make the RHS of Eq.(6) a source term, rather 
than a flux term, in order to apply the inlet flux boundary condition properly. The above transport PDE, Eq.(6), 
represents Fickian diffusion 6 in conservative form and assumes that all mechanical dispersion effects are lumped 
together with molecular diffusion in the axial dispersion term. 

The fifth, sixth, and seventh PDEs are coupled and solve for the sorbent temperature, T s , the gas temperature, T g , 
and the wall housing temperature, T c , respectively. Together these are referred to as the ‘Thermal Balance’ 
equations. The General Form PDEs are: 


dTg_ 

6 sorbentPg c pg + q x 


Kgxt sorbent q x j ™ l sg\ l s 1 

dT c d / dT c \ p,h gc (j g ■ 

PcCpc dt dx\ c dx) A c 


Again, all of the terms on the RHS of Eqs. (10-12) are source terms in COMSOL. Since u s and u w are identically 
zero everywhere (thus, u=u g explicitly throughout this document), no explicit boundary conditions are required for 
Eqs. (10 and 12). In Eqs. (10-12), subscripts s, g, A, and c refer to properties of the sorbent, gas mixture, ambient 
environment, and can housing, respectively. P/ = 27iR B , Po= 27 r(R B +d), and A = 7i((R B +d) 2 -R B 2 ) are the can inner 
perimeter, outer perimeter, and cross sectional area, respectively. Heat capacity, thermal conductivity, and thermal 
transfer coefficients are given by c p , k , and h, respectively. The heat of adsorption for a given sorbate/sorbent pair is 
given by dH and is assumed constant (see Table 1 below). The ratio of sorbent area to volume assumes spherical 
pellets and is given by A = i\-e sorhent )6AjJD veUets . The outlet boundary condition for Eq.(ll) is given by a zero 
gradient Constraint: 


dj, 

dx 


X=B 


= 0(13), 


while the inlet is given by a flux boundary condition: 


dTg 
gx dx 


X = 0 


~ Pg c pg u {^gin(t) ^ g ) ( 14 ). 
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The initial conditions for the c and q PDEs correspond to equilibrium loading for a small initial sorbate partial 
pressure, typically 1-5 Pa. The initial velocity is u in everywhere, while the initial pressure is set to P in -AP*x/Z?. All 
temperatures are initially uniform and set to the measured (but offset) T^ n (£ = 0). 

The boundary conditions described above for P, T, and c are all the measured values from a given test, with T 
and c being time-dependent. This is necessary since, as will be seen below, the variation from test to test (or even 
within a single test), even using the same packed bed and nominal test conditions for flow rate and inlet sorbate 
partial pressure, is often larger than the uncertainty in the model. 


1 . Dimensionless Numbers 

Many of the physical parameters used in Eqs. (10-12) are determined from correlations based on dimensionless 
quantities. The empirical relationships used here are appropriate for 1-D models in the regime of the CBT. The 
Schmidt number, Sc, pellet Reynolds number, Re, Peclet number 7 , Pe, Prandtl number, Pr, gas to sorbent Nusselt 
number, Nu G $, and gas to can Nusselt number, Nu GC 9 , are given by: 


Sc 


Re = 


Pg^AB 
u PgDpellets 


( 15 ) 


(16) 




0.73e, 


sorbent 


ReSc 


+ - 


f(x,r ) (l "f 

P c pg 


Pr = 


13 * 0-73e SO rbent \ 
ReSc )/ 

(18) 


Nu gs = 2 + 1.1 Pr 1/3 Re 06 (19) 
Nu gc = 2.63 Re°- 8 e~ 3D P ellets/R B (20). 


(17) 


In Eq.(17), f( x>r ) is 2 for axial ( Pe x ) and 10 for radial (. Pe r ) parameters. In Eq.(15), D AB is the binary mass diffusion 
coefficient 10 for either water vapor or carbon dioxide in N 2 as a function of P and T g . 

2. Toth 

The loading equation uses the Toth isotherm relationships: 


q* 


aP 

u-i VCi p 

b = b 0 e E/T Q (22) 
a = a 0 e E / T 3 (23) 

t T = to+^ (24), 

1 n 


( 21 ) 


where P vap = c R g T g is the sorbate partial pressure and a 0 , b 0 , t 0 , c 0 , and E are Toth coefficients for a given 
sorbent/sorbate pair. The coefficients used 12 are listed in Table 1, together with the heat of adsorption 12,13 for each 
pair. It is known that dH depends on loading, decreasing significantly as equilibrium loading is approached; the 
quantitative dependence, however, is not well known for many sorbents, so the impact is left for future work. For 
the silica gel sorbents, Grace Grade 40 and Sylobead B125, the same Toth parameters are used, even though they are 
experimentally derived from the former. The same is true for the 5 A zeolite sorbents, Grace Grade 522 and RK-38. 
There are no CBT experiments as yet with H 2 0/5A or C0 2 /SG systems. 
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Table 1. Adsorption Parameters for Sorbent/Sorbate Pairs 


Sorbate/Sorbent 

System 

a 0 

mol kg" 1 kPa" 1 

b % 

kPa' 1 

E 

K 

to 

Co 

K 

dW 

kJ mol" 1 

C0 2 /5A 

9.875xl0" 7 

6.761xl0" 8 

5.625xl0 3 

2.700xl0 _1 

-2.002x10' 

-38.0 

H 2 0/5A 

1.106xl0" 8 

4.714xlO" 10 

9.955xl0 3 

3.548xl0 _1 

-5.114x10' 

-45.0 

C0 2 /13X 

6.509xl0' 3 

4.884xl0' 4 

2.991xl0 3 

7.487xl0' 2 

3.805x10' 

-40.0 

H 2 0/13X 

3. 634x1 O' 6 

2.408x1 0" 7 

6.852xl0 3 

3.974x10"' 

-4.199 

-55.0 

co 2 /sg 

7. 678x1 O' 6 

5.164xl0" 7 

2.330xl0 3 

-3.053x10"' 

2.386xl0 2 

-40.0 

h 2 o/sg 

1.767xl0 2 

2.787xl0' 5 

1.093xl0 3 

-1.190xl0" 3 

2.213x10' 

-50.2 


3. Material Properties 

Here the expressions for material properties are given. Some are from COMSOL’s material libraries and some 
use the correlations in the previous section. Temperatures in the following expressions are in Kelvin and the derived 
values are in MKS units. 

The axial mass dispersion coefficient used in Eq.(6) for the transport of the sorbate through the bed is given by: 

D x = MPpe “ ets (25). 

P e x 

For the gas mixture, COMSOL expressions for N 2 are used: 

k g = 3.6969697xl0” 4 + 9.74353924x10"% - 4.07587413xl0"% 2 + 7.68453768xl0" 12 7/ (26) 

kgx = c pgPgD x ( 27 ) 

c pg = 1088.22121 - 0.3659419197^ + 7.88715035xl0"% 2 - 3.749223xl0"% 3 
+ 3.17599068xl0" 11 7/ (28) 

p. = 1.77230303xl0" 6 + 6.27427545x10"% - 3.47278555xl0" 1 % 2 + 1.01243201xl0" 1 % 3 (29). 

The effective axial gas thermal conductivity, k gx , comes from a similarity assumption 6 , which may not be valid in 
all regimes. The resulting values are typically as much as 50% larger than when using other, more complicated 
expressions 14 .; the impact of different k gx expressions on the breakthrough curves discussed here is minor. 

For the can properties, COMSOF functions for A1 6061 are used: 

k c = 30.71577 + 0.87509487; - 0.0021423317/ + 2.05104xl0"% 3 ,T C < 300 
= 79.478 + 0.33550017; - 2.701638xl0" 4 7/, T c > 300(30) 
c pc = -5.605966 + 6.7447257; - 0.015359127/ + 1.619759xl0"% 3 - 6.214813xl0"% 4 (31) 
p c = -2736.099 - 0.009406247^ - 6.040342xl0" 4 7; 2 + 8.988964xl0" 7 7; 3 - 5.405225xl0" 1 °7’ 4 (32) 

Note the piecewise expression used for k c . 

For the sorbent properties, since they are not well known, constant values are used 12,15,16 . They are listed for the 
sorbents used in this work in Table 2. For Sylobead B152, the sorbent density and heat capacity are assumed to be 
the same as for Grace Grade 40. Also listed in Table 2 are the mean pellet diameter, total dry sorbent mass, and 
resulting mean porosity for the sorbent packings used in the CBT experiments. The heat capacity of the pellets can 
change significantly as they get loaded. Thus, using, from COMSOF, the heat capacities of water vapor: 

c pH20 = 1.244444465x10 " 10 T 5 6 - 2.1702974711043xl0" 7 T s 5 + 0.000157628221317181 T s 4 
- 0.0610285680460429 T s 3 + 13.28450401510987/ - 1541.532490628437; 

+ 74503.5945678409 (33), 


And of carbon dioxide: 

c pC02 = 459.9 1 3 2 5 8 + 1.86 4 8 7 9 9 67; - 0.002129215197/ + 1.22488004xl0" 6 7/ (34), 
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a mass-weighted loading-dependent effective heat capacity can then be used for c ps : 

_ ( C p(H20,C02) C lM( < H20 > C02)+ C psoPs) 

Cps " (gM (H20 , C02) + Ps ) ^ 

Eq.(35) typically reduces the pellet heat capacity by -10% when the sorbate is H 2 0, but is insignificant for C0 2 . 
The thermal transfer coefficients for the gas are determined from Nu: 


hsg ~ NUcsPg/Dpeiiets ( 36 ) 


hgc 


NU GC Pg 

(2 R b ) 


(37), 


while for the ambient to can heat transfer, a nominal constant small value of h Ac - 0.1 W/m 2 /K was used, since the 
insulation was not firmly bonded with the cylinder. As a result, the value used for T A , the ambient temperature, is 
not very critical, so a constant 19°C value, representative of the CBT laboratory, was used in this work. However, 
varying T A to the measured value for each test (when available) would result in slightly improved late time 
temperature profiles (see below). Also, h sg is typically -150 W/m 2 /K, so the gas and sorbent are well thermally 
coupled, while h gc is typically -25 W/m 2 /K, so the can is less coupled to the flow. As a result, in the results 
discussed below, typically, I T g - T s \ «1°C and IT^-T^I »1°C. 


Table 2. Properties of Sorbents as Packed in the CBT 


Sorbent 

Ps 3 

kg m 

c ps0 (@26°C) 
J kg' 1 K' 1 

k s 

W m 1 K 1 

Dpellets 

mm 

mass 

g 

£ sorbent 

5A (Grace Grade 522) 

1190 

750 

0.152 

2.22 

125.0 

0.331 

5A (RK38) 

1370 

650 

0.144 

2.10 

119.3 

0.445 

13X (Grace Grade 544) 

1260 

800 

0.147 

2.19 

107.4 

0.457 

SG (Grace Grade 40) 

1240 

870 

0.165 

2.90 

111.7 

0.415 

SG (Sylobead B125) 

1240 

870 

0.151 

2.25 

127.0 

0.348 


C. COMSOL Results 

The only free parameter required to fit the CBT data is k m , the mass transfer coefficient used in the LDF model 
(see Eq.(7)). It is not expected to be sensitive to different test conditions, such as flow rates, vapor pressure, or 
temperatures. Ideally, it is only a function of the sorbent/sorbate pairing, but in practice it can vary due to, for 
example, geometry differences when R B is only a few times D pe p ets , making arbitrary predictive applications with the 
LDF model problematic. Here, the same value of k m is used for all CBT experiments of a given sorbent/sorbate 
pairing. Note that k m drives the slope of the vapor pressure rise curve, so for the model to be self-consistent, it 
should reproduce that slope; that is, even k m is not a true free parameter even for a single test of a given 
sorbent/sorbate pairing. The intent is to derive k m for a given sorbent/sorbate pair and use it in future predictive 
modeling of CDRA-related systems. It is to be emphasized that once k m is determined for a given test of a 
sorbate/sorbent system, the other tests of that sytem (at different flowrates and partial pressures) are predictively 
modeled. 

In all tests, there was a subsampler used to measure water dew point or carbon dioxide partial pressure; the 
volumetric flow rate for the subsamplers was fixed at 0.7 LPM. As a result, the upstream flow was set to the 
intended flow rate plus 0.7 SLPM, so the flow entering the sorbent bed should be within -10% of the intended 
SLPM. Due to using standardized flow meters, the actual volumetric flowrate into the bed is typically -10% higher 
than the SLPM values listed for these experiments. 

The COMSOL results for the predicted and experimentally measured exit temperatures and vapor pressures are 
shown in the following figures. In some cases, multiple nominally identical tests were run, while in others only a 
single test is discussed. Sometimes multiple flowrates and inlet partial pressures were tested, sometimes only one 
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condition was tested. The derived values of k m used for the various sorbent/sorbate systems are listed in Table 3. It 
can be seen that k m varies by a factor of 5, with 0.002 s' 1 being a typical value. 

Table 3. COMSOL Derived Mass Transfer Coefficients 


Sorbent/Sorbate 

k m 

System 

s" 1 

Sylobead/H 2 0 

0.002 

Grade 40/H 2 O 

0.00125 

Grade 544/H 2 0 

0.0007 

RK38/C0 2 

0.003 

Grade 522/CQ 2 

0.0035 


H 2 0 on Silica Gel Grade 40 (F=8SLPM 0.5C Dew Point) 
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H 2 0 on Silica Gel Grade 40 (F=8SLPM 10C Dew Point) 



H 2 0 on Silica Gel Grade 40 (F=16SLPM 0.5C Dew Point) 
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H 2 0 on Silica Gel Grade 40 (F=16SLPM 10C Dew Point) 



H z O on Silica Gel Grade 40 (F=8SLPM 0.5C Dew Point) 



time (min) 

Figure 7. Vapor pressures for SG Grade 40 8 SLPM 0.5°C Dew Point. 
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H 2 0 on Silica Gel Grade 40 (F=8SLPM 10C Dew Point) 
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Figure 8. Vapor pressures for SG Grade 40 8 SLPM 10°C Dew Point. 
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Figure 9. Vapor pressures for SG Grade 40 16 SLPM 0.5°C Dew Point. 
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H 2 0 on Silica Gel Grade 40 (F=16SLPM 10C Dew Point) 



Figure 10. Vapor pressures for SG Grade 40 16 SLPM 10°C Dew Point. 


H 2 0 on Silica Sylobead B125 (F=16SLPM 10C Dew Point) 
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H £ 0 on Silica Sylobead B125 (F=16SLPM IOC Dew Point) 

7 

6 

5 

4 

3 

2 

1 

0 

0 200 400 600 800 1000 1200 1400 

time (min) 

Figure 12. Vapor pressures for SG Sylobead 16 SLPM 10°C Dew Point. 



H 2 0 on Zeolite 13X Grade 544 (F=9.5SLPM 10C Dew Point) 



time (min) 

Figure 13. Temperatures for 13X Grade 544 9.5 SLPM 10°C Dew Point. 
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H 2 0 on Zeolite 13X Grade 544 (F=9.5SLPM IOC Dew Point) 



1. H 2 0 on Silica Gel Grade 40 

Two flowrates, 16 SLPM and 8 SLPM, and two inlet partial pressures, corresponding to dew points of 0.5°C and 
10°C, were run with silica gel Grace Grade 40 sorbent. The two flowrates correspond to Reynolds numbers of ~ 150 
and ~70, respectively. At least 2 tests were run at each condition, though only two tests for each condition are 
discussed here. Figures 3 to 6 show the temperatures and Figs. 7 to 10 show the vapor pressures. The results show 
significant variation between and within tests (see, e.g., Fig. 3 and Fig. 7). The most significant temperature 
difference between the COMSOL models and the experiments appears as a ‘hump’ at low flowrate and low dew 
point (Fig. 3). The hump is visible in Fig. 5 as well, but is much smaller in amplitude. The hump is due to too large 
of an initial rise in temperature (also evident in Fig. 4), followed by a near-leveling off that is not seen in the data, 
though Fig. 5 only shows the latter. All of this suggests issues with the correlations breaking down at low flowrate 
and/or vapor pressures. The partial pressure results show that the data at late times do not approach the inlet values; 
the reason for this consistent deficiency is unknown but may be related to calibration issues with the dew point 
sensors, since it does not occur in the C0 2 tests (see below). Also, some of the temporal variation in the data are 
real (reflected in simultaneous changes in temperature and vapor pressure), while some are not; in the models, all are 
taken as real so that some variations (see, e.g., Fig. 8) are incorrectly propagated (but for the right reasons). The 
modeled temperature rise 1 inch inside the bed inlet is slower than observed in the experiments (not shown in the 
figures) due to the flux boundary condition resulting in the inflowing sorbate concentration not immediately rising to 
the experimental inlet value; this is true in all of the H 2 0 CBT models to various degrees. 

2. H 2 0 on Silica Gel Sylobead B125 

Only one flow condition was tested for the Sylobead B125 sorbent: 16 SLPM and 10°C dew point. Figs. 11 and 
12 show the model and experiment results for the temperature and water partial pressure, respectively. The 
comparison is not as good as for Grade 40. This is likely due to two factors. First, the Sylobeads have a 
significantly different size and shape distribution, such that the packed mean porosity is significantly different (see 
Table 2) and thus resulting in a different flow field in the bed. Second, the used density, heat capacity, and Toth 
parameters were derived for Grace Grade 40, not the Sylobead B 125 sorbent. 
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C0 2 on 5A Grade 522 (F=8 SLPM; 5 Torr partial pressure) 



Figure 15. Temperatures for 5A Grade 522 8SLPM 5 Torr partial pressure of C0 2 . 


C0 2 on 5A Grade 522 (F=8 SLPM; 2.5 Torr partial pressure) 
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C0 2 on 5A Grade 522 (F=8 SLPM; 5 Torr partial pressure) 


CC 

CL 


> 

CL 



CO z on 5A Grade 522 (F=8 SLPM; 2.5 Torr partial pressure) 
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C0 2 on 5A Grade 522 (F=16 SLPM; 5 Torr partial pressure) 



C0 2 on 5A Grade 522 (F=16 SLPM; 2.5 Torr partial pressure) 
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C0 2 on 5A Grade 522 (F=16 SLPM; 5 Torr partial pressure) 



C0 2 on 5A Grade 522 (F=16 SLPM; 2.5 Torr partial pressure) 
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3. H 2 0 on 13X Zeolite Grade 544 

Only one CBT experiment was conducted with H 2 0 on Grace Grade 544 due to difficulties in regenerating the bed. 
The inlet pressure, pressure drop, and ambient conditions were not recorded; they are assumed to be comparable to 
the other tests. The inlet water vapor partial pressure was 1.2 kPa, corresponding to a dew point of 10°C. Although 
the flowrate was supposed to be 16 SLPM, the flow meter reference settings and calibrations may have been off for 
this test. Thus, the results given here in Figs. 13 and 14 assume an actual flowrate of 9.5 SLPM. This gives a good 
fit to the water vapor rise curve (Fig. 14). The temperature peak in the model is too low and drops off too slowly 
compared to the data. A combination of higher heat of adsorption and poor insulation would be required to explain 
this. However, another explanation, instead of a low flowrate, is that the porosity was substantially incorrect, due to, 
for example, in incorrect bulk pellet density or gross non -sphericities in the pellets, as well as the pellets starting 
with an initially extremely low loading (as if coming straight from bake out with no N 2 flow to equilibrate the 
temperature in the bed). 

4. CO 2 on 5 A Zeolite Grace Grade 522 

It is well known that C0 2 competes with H 2 0 on zeolites for adsorption, with water being much more easily 
adsorbed. Thus, for the C0 2 CBT experiments, since the purge gas is not perfectly dry (the purge has a dew point of 
— 70°F), a larger initial loading is assumed in order to mimic the H 2 0 loading that would occur during the N 2 purge. 
For Grade 522, an initial uniform loading of 350 mol/m 3 was thus assumed. This is an approximation based on the 
equilibrium water loading of Grade 522, but, as a result, the model does not properly capture the initial partial 
pressure breakthrough curve. Two flowrates, 16 and 8 SLPM, and two inlet partial pressures, 2.5 and 5.0 Torr, were 
run for C0 2 on 5 A Grade 522 sorbent. Multiple tests were run at each condition except for 16 SLPM and 5 Torr. 
Comparisons between the experiments and the COMSOL models are shown in Figs. 15 through 22. As expected, 
the model consistently misses the gentle initial rise of the ‘heel’ of the C0 2 partial pressure curves, particularly at 
low partial pressures. In addition to the artificial initial loading, this is also likely due to the 1-D model not being 
able to capture the channel flow near the walls, where breakthrough will occur more quickly; C0 2 adsorption is 
more sensitive to this effect than H 2 0 since it occurs more than ~5x more quickly (compare, for example, Figs. 4 
and 15). A better approach may be to turn off loading in some initial depth of the bed where water prevents any 
C0 2 loading; this is left for future work since the precise details of the purge (duration, sorbate partial pressure, 
volumetric flowrate) are not known for the purge gas used in the CBT experiments. In addition, the model peak 
temperatures, particularly for the high flow rate tests, are consistently too low. This may be due to an inaccurate 
heat of adsorption or the fact that the Nu correlations used to estimate thermal coefficients were based primarily on 
H 2 0 data, not C0 2 data. For the C0 2 tests on both Grade 522 and RK38 (see below), the sorbent pellets were sieved 
with an 8x12 sieve, reducing the spread of particle size in the bed to between 2.38 and 1.41mm. Other tests using a 
20x50 sieve (not discussed here further) for Grade 522 showed very similar results. 

5. CO 2 on 5 A Zeolite RK38 

The initial loading for the C0 2 on RK38 tests was set to 250 mol/m 3 . Two flowrates, 16 and 8 SLPM, and two 
inlet partial pressures, 2.5 and 5.0 Torr, were run. Multiple tests were run at each condition except for 16 SLPM and 
2.5 Torr. Comparisons between the experiments and the COMSOL models are shown in Figs. 23 through 30. 
Again, as expected, the model consistently misses the rise of the C0 2 partial pressure curves. In addition, due to the 
initial non-equilibrium loading, there are some numerical instabilities in the sorbate concentration at initial times 
which show up as wiggles in the 5 Torr partial pressure plots (see, e.g., Fig. 25). Also, as for Grade 522, the model 
peak temperatures are too low, particularly at high flow rates. This suggests that issues with the C0 2 model or 
testing procedure in general are more likely the problem than sorbent-specific properties such as the heat of 
adsorption. 

The adsorption capacity, the moles of sorbate per unit mass of sorbent that the system can adsorb, was measured 
for the 16 SLPM 5 Torr flow condition tests to be 1.15+0.05 mol/kg. In the models, the loading throughout the bed 
is within a few percent of the equilibrium loading by the end of the tests. Thus, the calculated theoretical capacity is 
qd p s ; it falls within the uncertainty range of the measured values. This argues against any significant fraction of the 
sorbent in the bed being ‘inactive’ due to the presence of water. Another possibility is that all of the C0 2 Toth 
isotherms for zeolites 11 were derived with sorbents that were contaminated with some unknown amount of H 2 0 due 
to too low of an activation temperature. This would make them inaccurate at low C0 2 loading levels. 


21 

International Conference on Environmental Systems 



22 


C0 2 on 5A RK38 (F=8 SLPM; 5 Tore partial pressure) 



time (min) 

Figure 23. Temperatures for 5A RK38 8 SLPM 5 Torr partial pressure of C0 2 . 


C0 2 on 5A RK38 (F=8 SLPM; 2.5 Torr partial pressure) 



time (min) 

Figure 24. Temperatures for 5A RK38 8 SLPM 2.5 Torr partial pressure of C0 2 . 
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C0 2 on 5A RK38 (F=8 SLPM; 5 Torr partial pressure) 



C0 2 on 5A RK38 (F=8 SLPM; 2.5 Torr partial pressure) 



time (min) 

Figure 26. Partial pressures for 5A RK38 8 SLPM 2.5 Torr partial pressure of C0 2 . 
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C0 2 on 5A RK38 (F=16 SLPM; 5 Torr partial pressure) 



Figure 27. Temperatures for 5A RK38 16 SLPM 5 Torr partial pressure of C0 2 . 


C0 2 on 5A RK38 (F=16 SLPM; 2.5 Torr partial pressure) 



Figure 28. Temperatures for 5A RK38 16 SLPM 5 Torr partial pressure of C0 2 . 
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C0 2 on 5A RK38 (F=16 SLPM; 5 Torr partial pressure) 



C0 2 on 5A RK38 (F=16 SLPM; 2.5 Torr partial pressure) 



time (min) 

Figure 30. Partial pressures for 5A RK38 16 SLPM 2.5 Torr partial pressure of C0 2 . 
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IV. Conclusion 

Using COMSOL, we have been able to derive a predictive model for adsorption physics in a variety of flowrate 
and sorbate partial pressure regimes and for a variety of sorbents. For water vapor on silica gel and zeolites, the 
limitations are in general due to the experimental data, not the model, due to variations and uncertainties in the dew 
point measurements, the ambient conditions, and unsteady flow conditions. For carbon dioxide on zeolites, the 
model does not compare as well to the data. However, it is not clear if this is a deficiency of the model, such as the 
LDF, which uses a single constant parameter to determine the loading rate, or of the experiments, due to competing 
pellet loading from water vapor in the purge gas. Future C0 2 tests will need to carefully track the amount of H 2 0 
loading on the sorbent for any given test. In addition, the 1-D nature of the model is not capturing some of the 
physics due to channeling of the carrier gas near the walls of the sorbent bed. Future work will extend the PDEs 
used here to 2-D axisymmetry and will include the impact of using a radially dependent porosity due to packing. 
Since no purge gas is perfectly dry, another extension of this work is to include binary Toth relationships, 
particularly for the C0 2 , in order to properly capture the competition between H 2 0 and C0 2 . Also, it is known 13 that 
the heat of adsorption decreases with increasing loading; this effect, when known, will be incorporated in future 
models. The values of k m found in this work will be used in broad applications of sorbent/sorbate systems to 
estimate performance without the need of testing. 
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